#include <iostream>
using namespace std;
#include "diff_equ_1D_gas_phase_rec.h"

int main(int argc, char** argv)
{
	diff_equ_1D_gas_phase_rec	GasPhase(0.5,0.4,133.32,300,850,0.01,0.00001,10,0.005,0.01,0.1,"test_Xin01.dat",1);
	for(int j=0;j<1;j++)
	{
		GasPhase.reset();
		cout << "Run #" << j+1 << "\n";
		while(GasPhase.notfinish())
		{
			GasPhase.update(0.001);
			cout << GasPhase.time << "\n";
		}
		cout << "\n\n\n";
	}
	cout << "Test DIFFUSION_equ_1D success!\n";
	system("pause");
}